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Abstract 

Coronagraphic imagery of the circumstellar disk around HD 169142 in //-band 
polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered 
by dust particles at the disk surface in 0."2 < r < l."2, or 29 < r < 174 AU, is successfully 
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detected. The azimuthally-averaged radial profile of the PI shows a double power-law 
distribution, in which the Pis in r = 29 — 52 AU and r = 81.2 — 145 AU respectively 
show r^ 3 - dependence. These two power-law regions are connected smoothly with a 
transition zone (TZ), exhibiting an apparent gap in r = 40 — 70 AU. The PI in the 
inner power-law region shows a deep minimum whose location seems to coincide with 
the point source at A = 7 mm. This can be regarded as another sign of a protoplanet 
in TZ. The observed radial profile of the PI is reproduced by a minimally flaring 
disk with an irregular surface density distribution or with an irregular temperature 
distribution or with the combination of both. The depletion factor of surface density 
in the inner power-law region (r < 50 AU) is derived to be > 0.16 from a simple model 
calculation. The obtained PI image also shows small scale asymmetries in the outer 
power-law region. Possible origins for these asymmetries include corrugation of the 
scattering surface in the outer region, and shadowing effect by a puffed up structure 
in the inner power-law region. 

Key words: stars: individual (HD 169142) — stars: pre- main- sequence — plan- 
etary systems: protoplanetary disks — planetary systems: planet-disk interactions 
— infrared: planetary systems — infrared: stars 


1. Introduction 

Circumstellar disks around young stars are the key targets to obtain a better under- 
standing of the formation of planetary systems (Dullemond & Monnier 2010; Williams & Cieza 
2011). During ongoing formation of a planet, the disk is expected to manifest observable sig- 
natures inside, such as an inner hole, a gap, spirals and other asymmetric patterns. Recent 
progress in the technique of differential imaging at near infrared wavelengths finally allows us 
to detect small-scale complex features in the disks (e.g., Hashimoto et al. 2011; Hashimoto et 
al. 2012; Muto et ah 2012; Mayama et ah 2012; Grady et ah 2013; Follette et ah 2013). Most 
of these objects are classified as transitional disks, which show significant amount of excess 
emission at wavelengths A > 10 /irn but little excess at shorter wavelengths, suggesting that 
dust in the inner regions has already been depleted (Strom et ah 1989; Skrutskie et ah 1990). 
Imaging surveys of transitional disks in dust thermal emission at submillimeter wavelengths 
clearly show that these disks do commonly possess inner holes (Brown et ah 2009; Andrews 
et ah 2011). Recently, even more complex distributions of dust emission on a smaller size 

* E-mail: momose@mx.ibaraki.ac.jp 
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scale are also revealed in some of these objects (e.g., Casassus et al. 2013; van der Marel et 
al. 2013; Fukagawa et al. 2013; Perez et al. 2014). These asymmetric features may be caused 
by trapping of dust particles in a vortex with a higher pressure, which can lead to efficient 
formation of rocky planetesimals (Regaly et al. 2012; Birnstiel, Dullemond & Pinilla 2013). 

HD 169142 is a Herbig Ae star located at 145 pc from the Sun (Sylvester et al. 1996), 
and the spectral type of the star is classified as A5Ve (Dunkin, Barlow & Ryan 1997) or A8Ve 
(Grady et al. 2007). The age is estimated to be 61® Myr from the ages of two pre-main sequence 
stars located 973 to the southwest and are comoving with HD 169142 (Grady et al. 2007). The 
luminosity of 9.4L 0 is used by Meeus et al. (2012), while Manoj et al. (2006) derive 18.2L 0 
with Ay = 0.61 (see also Osorio et al. 2014). Judging from the observed B — V color 0.26, 
however, Ay fa 0.5 was the case only if the stellar spectral type was as early as B8. Wagner 
et al. (2015), in fact, find a good fit to the SED, including the UV data, for the spectral type 
of A7, the effective temperature of 7500K and the stellar radius of 1.7 Rq (resulting in 8.7 L 0 ) 
with no foreground extinction. It exhibits significant amount of excess emission at longer 
wavelengths and its spectral energy distribution (SED) is categorized as group I (Meeus et al. 
2001). Interferometric observations of CO lines at millimeter wavelengths revealed a gas disk 
with the radius of 2" (fa 300AU), and their velocity fields are consistent with Keplerian rotation 
of the central stellar mass 2 M 0 and its inclination angle 13° (Raman et al. 2006; Panic et al. 
2008). Scattered light at near infrared wavelengths from the disk outer regions was also imaged 
with coronagraph instruments (Grady et al. 2007; Fukagawa et al. 2010). 

More recently, the inner regions of the disk have been investigated by both the analysis 
on its broadband SED and higher resolution imagery. It was claimed that there should be an 
inner cavity around HD 169142 to reproduce its SED, but its size was rather unclear solely 
from the SED analysis (Grady et al. 2007; Meeus et al. 2010). Honda et al. (2012) carried out 
imaging of the disk at A = 18.8 and 24.5 //m with Subaru/COMICS, and by combining these 
data with a two-dimensional radiative transfer model for the SED, they determined the location 
of the inner wall of the disk to be 23 ± 3 AU from the star. They also speculate that the SED 
classification as group I may correspond to the object accompanied by a (pre)transitional disk 
in which there is a large inner hole and the inner wall of the disk just outside the hole can 
account for its characteristic large amount of excess emission at mid-infrared wavelengths (see 
also Maaskant et al. 2013). Quanz et al. (2013) made observations of polarization intensity 
(PI) in H - band (A = 1.65 /im) with VLT /NACO . They found a bright ring of r « 25 AU that 
agreed with the estimate from the SED analyses (Meeus et al. 2010; Honda et al. 2012), as well 
as an annular gap extending r = 40 — 70 AU. Furthermore, imaging of dust thermal radiation 
at A = 7 mm with J-VLA also shows asymmetric ring-like emission at r fa 25 AU and a hint 
of an outer gap at 40 — 70 AU, suggesting that the structure seen in the scattered light at 
near infrared wavelengths may correspond to real variations in surface density (Osorio et al. 
2014). Although these are not the decisive evidence for the existence of protoplanet (s), such a 
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gap is reminiscent of that induced by a planet in the disk (e.g., Lin & Papaloizou 1993; Crida, 
Morbidclli & Masset 2006). The latest observations in L'-band (A = 3.8 pm) with an annular 
groove phase mask (AGPM) vector-vortex coronagraph mounted on VLT detect a companion 
candidate inside the inner cavity (07156, or 23 AU from the star), which might be the direct 
detection of radiation from a protoplanet (Biller et al. 2014; Reggiani et ah 2014). 

This paper presents coronagraphic imagery of the circumstellar disk around HD 169142 
in H - band polarized intensity (PI) with Subaru/HiCIAO. We describe the details of the ob- 
servations in §2 and present the results in §3. Based on our results as well as those obtained 
by previous studies, we discuss in §4 the disk structure and estimate the nature of a possible 
protoplanet embedded in the disk from our results. 

2. Observations and Data Reduction 

2.1. Observations with Subaru/HiCIAO 

HD 169142 was observed in //-band (A = 1.65 pm) with the high-contrast imaging 
instrument HiCIAO (Tamura et ah 2006; Hodapp et ah 2008; Suzuki et ah 2010) on the Subaru 
Telescope on 2011 May 23 UT as part of Strategic Explorations of Exoplanets and Disks with 
Subaru (SEEDS; Tamura 2009). The data were taken in a combined angular differential imaging 
(ADI) and polarization differential imaging (PDI) mode with a held of view of 10" x 20" and a 
pixel scale of 9.5 milli-arcsecond (mas). A circular occulting mask 0715 in radius was used to 
suppress the bright stellar halo. The half-wave plates were periodically placed at four angular 
positions from 0°, 45°, 22.5° and 67.5° in sequence with 30s exposure per wave plate position. 
The total integration time of the PI image was 1080s after removing low quality frames with 
large FWHMs by careful inspection of the stellar point-spread function (PSF). Data of the 
reference star HD 166903, whose apparent magnitude in //-band is 5.837 mag (SIMBAD), were 
also taken in the same observing mode. It was proved that a stable stellar PSF of FWHM = 
0706 was achieved for the selected images of HD 169142 by the adaptive optics system A0188 
(Hayano et ah 2004; Minowa et ah 2010). 

2.2. Data reduction 

The data reduction procedure was the same as that described by Muto et al. (2012). 
The raw images were first corrected with IRAF 1 for bias and dark current, and hat-fielding 
was performed after sky subtraction. We applied a distortion correction using globular cluster 
M5 data taken within a few days, by IRAF packages GEOMAP and GEOTRAN. The images 
of Stokes (Q,U) parameters were obtained in the standard way (e.g., Hinkley et al. 2009); By 
subtracting two images of extraordinary and ordinary rays at each wave plate position, we got 

1 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 
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+Q,—Q,+U, and —U images, from which 2 Q and 2 U images were made by another subtraction 
to eliminate remaining aberration. Instrumental polarization of HiCIAO was corrected by the 
way described by Joos et al. (2008). 

We needed to fix the rotation center of the field of view by ADI in de-rotating the Q and 
U images. It was determined in the following way. A pair of /-frames that were taken more 
than 300s apart from each other was picked up, and then the rotation center and rotation angle 
for the pair were derived so that the two background stars in the differential image between 
the two frames were best cancelled out. These values were collected for all the possible /-frame 
pairs, and their average position was regarded as the image rotation center in our observations. 
There were four background stars in each frame, but two of them were not used because their 
positions were so far from the image center that the distortion correction might be imperfect. 
It proved that the mask center position was shifted from the image rotation center by (+1,-4) 
pixels in (a, h), or (+9.5 mas, —38 mas). On the other hand, the stellar position was estimated 
from elliptical fitting to an iso-intensity contour of the periphery of stellar halo in /-frames. 
The stellar location measured in this way was proved to be the same as the image rotation 
center. The Q and U images were de-rotated at this center, then averaged to increase the 
signal-to-noise. PI was then obtained by PI = sJQ 2 + U 2 . The misalignment between the mask 
center and image rotation center made the resultant working angle > 0/2 from the star, which 
is slightly larger than the actual mask size. 

3. Results 

Figure 1 shows the PI image of HD 169142 in H - band with Subaru/HiCIAO. The infrared 
emission scattered by dust particles at the disk surface is successfully detected. The typical PI 
value is ~ 0.21 mJy arcsec” 2 in r > 3", well beyond the disk outer radius derived from the CO 
observations (Raman et al. 2006; Panic et al. 2008) and the scattered light of the disk should be 
negligible. This corresponds to the la level of 0.168 mJy arcsec” 2 if the values in the region of 
no polarized emission follow the Rayleigh distribution. In the following quantitative analyses, 
we will regard the PI emission of > 0.5 mJy arcsec” 2 (= 3<r level) as a significant detection. The 
region of significant detection extends up to r ~ 1/2 from the star, as presented in figure lc. 
Figure Id shows the polarization vectors derived from Stokes Q and U images. The polarization 
vectors align azimuthally, as expected for the stellar light scattered at the disk surface. 

The overall structure of the PI is roughly consistent with that obtained by Quanz et al. 
(2013) with VLT/NACO . The obtained PI in figure la shows an axisymmetric distribution 
as follows: (i) the bright ring just outside the mask, (ii) a ring-like feature at r « 0/55, or 
80 AU from the star, (iii) the apparent gap between these two ring-like features, and (iv) the 
outermost part in which the PI gets gradually weaker at a larger radius. Note that the PI is 
significantly detected in the gap (iii) with at least 9cr. These axisymmetric features can be seen 
even more clearly in the smoothed image convolved with a Gaussian function (figure lb). The 
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Fig. 1 . (a) Coronagraphic image of the polarization intensity (PI) in ff-band around HD 169142 obtained 
with Subaru/HiCIAO. White cross is the stellar position, and the green circle shows the regions affected 
by the mask, or <0.2" from the stellar position, (b) PI image convolved with a Gaussian function whose 
FWHM is 8 pixels, or 07072. (c) Same as (a), but in different color code so that asymmetric features 
are easily identified. White broken lines are the contours of 0.5 mJy arcsec -2 , corresponding to the 30- 
level (see the text). Solid lines in white and black are the contours at 1.0 mJy arcsec -2 (6 a) and 1.5 mJy 
arcsec -2 (9o-) levels, respectively. The contours are generated in the convolved image shown in (b). (d) 
Same as (b) , but in the color code as (c) , superposed on the polarization vectors where the PI is detected 
above the 3 <r level. The polarization vectors are drawn every 8 pixels. 
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inner bright rim at r ^ 072 (Quanz et al. 2013) is only marginally confirmed because the inner 
working angle of our observation is slightly larger than that of VLT/NACO. 

To describe the axisymmetric distribution of the PI more quantitatively, its azimuthally- 
averaged radial profile is presented in figure 2; the projection on the sky is corrected under the 
assumptions that the disk inclination angle is 13° and the disk major axis is in PA= 5° (Raman 
et al. 2006). The radial profile shows a double power-law distribution, as expressed by 

- 3.0025 
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29 AU 


- 3.0037 
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29 AU 

in 81.2 < r < 145 AU, respectively. The formal statistical uncertainties in the amplitudes 
are estimated to be 10% for (1) and < 1.3% for (2), and those in the exponents are <0.1 
dex. These two power-law regions are connected smoothly with a “transition zone” (hereafter 
denoted by TZ) extending r = 52 — 81 AU. The apparent gap located in r ^ 40 — 70 AU (Quanz 
et al. 2013; Osorio et al. 2014) consists of the outer part of the inner power-law region and TZ, 
while the bright ring at r & 0755 in the PI image (figures la and lb) corresponds to the inner 
boundary of the outer power-law region. It should be noted, however, that the surface brightness 
normalized by the r -3 power-law (figure 2c) changes smoothly even in TZ. There is no break 
of the exponent at r = 120 AU claimed by Quanz et al. (2013), though a subtle irregularity 
may be identified at r ^ 100 AU. The r -3 -dependence was also found in the total intensity 
distribution at 1.1 /irn in r > 80 AU with the coronagraphic imaging by HST /NICMOS (Grady 
et al. 2007). It is also consistent with the radial exponent in 116 < r < 174 AU derived from 
the imagery in 77-band with Subaru/CIAO (—3.0 ±0.2; Fukagawa et al. 2010), suggesting that 
there is no significant radial variation of polarization degree in these regions. The polarization 
degree estimated from the comparisons between the PI in figure 1 and the total intensity by 
Fukagawa et al. (2010) is ~ 65 % in 078 < r < 171. This is higher than that obtained in the 
77-band PI imaging of the disk around HD 142527 (Avenhaus et al. 2014). It may be due to the 
nearly pole-on configuration of HD 169142; the scattering angle at the disk surface should be 
~ 90° at which the polarization degree is ~ 100% when the particle size is much smaller than the 
wavelength (Bohren & Huffman 1983). Alternatively, it might be due to an underestimate of 
total scattered-light intensity caused by a systematic error during PSF subtraction by Fukagawa 
et al. (2010). The rest of the discussion in this paper will not be affected by any uncertainties 
in polarization degree. 

The PI image also manifests asymmetric features, which are easily identified in figure 
lc. To see the asymmetry at each radius, azimuthal distributions of the PI in annuli of 14.5 
AU width are presented in figure 3. In the innermost annulus (29.0 < r < 43.5 AU), the 
northwestern part (260° < PA < 340°) is quite bright. This may be a part of the inner rim 
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Radius [arcsec] 

0.2 0.4 0.6 0.8 1.0 



Fig. 2. (a) Averaged radial profile of the PI above the 3 a level (0.5 mJy arcsec -2 ), between (29—159.5) 
AU, or 0"2 — l."l. Averages and their uncertainties are calculated in annuli of 2.9AU (= 0"02) width after 
the correction for the projection on the sky under the assumptions that the position angle of the disk 
major axis is 5° and the disk inclination angle is 13°. (b) Results of the fitting by a power-law function 
to the data points in 29—52.2 AU (red) and in 81.2—145 AU (blue); the data points used in the power-law 
fitting are shown in blue and green dots. Radial profile in (a) is shown by gray dotted line, (c) Same as 
(a), but normalized by qjgAijj where r 2 gAU is a dimensionless radius normalized by 29 AU. 

distorted outwardly (Quanz et al. 2013). Except for this innermost annulus, all the azimuthal 
distributions have their maxima at PA^i 110°, close to the disk minor axis east to the star. This 
may imply that the eastern part is the near side of the disk and is brighter than the western 
side because of the forward scattering. In the outer power-law region, the PI near the major 
axis (PA = 5° and 185°) is dimmer than those near the minor axis (PA = 95° and 275°), which 
is also found in figure 4. As a consequence, the lowest contour of the significant detection in 
figure lc is slightly elongated in east-west direction. This situation is completely different from 
other disks in scattered light detected with Subaru/HiCIAO, in which the PI is brighter along 
the major axis (e.g., Muto et al. 2012; Kusakabe et al. 2012). The variation in each annulus 
is systematically larger at inner radii; In the inner power-law region and TZ, the amount of 
variation in each annulus is 4.5 — 3.3, but those of annuli in the outer power-law region is less 
than 2.6. 

To clearly identify these asymmetric features in the image, figure 5 shows the map of 
“modulated PI”, the PI normalized by r^Am where r 2 gAu is a dimensionless radius normalized 
by 29 AU. Two distinct power-law regions, divided by TZ where the contours are crowded, can 
be recognized. As indicated by the normalized intensities in equations (1) and (2), the average 
of modulated PI in the inner and outer power-law regions should be ~ 21 mJy arcsec -2 and 
~ 96 mJy arcsec -2 (see also figure 2c). However, the variation of modulated PI in the inner 
power- law region is quite remarkable, ranging from 10 to 40 mJy arcsec -2 . This is significant 
compared to the uncertainties estimated from the original PI map (3cr = 0.5 mJy arcsec -2 in the 
original PI corresponds to ~ 7.8 mJy arcsec -2 in modulated PI at r = 0."5). These asymmetries 
in the inner power-law region could be related to the perturbation by a protoplanet inside the 
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(a) Inner power-law region 


(b) Transition Zone 


(c) Outer power-law region 
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Fig. 3. Azimuthal distributions of the PI in annuli of 14.5AU (or 0"1) width. The projection on the sky 
is corrected under the same assumptions as those described in the caption of figure 2. The abscissa starts 
from PA= 5°, corresponding to the major axis of the disk, (a) in the inner power-law region, (b) in the 
transition zone (TZ), and (c) in the outer power-law region. 
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Fig. 4. (a) Radial profiles of the PI along the major axis of the disk: PA = 5° ± 5° in red and PA = 
185° ± 5° in blue, (b) Radial profiles of the PI along the minor axis of the disk: PA = 95° ± 5° in red and 
PA = 275° ± 5° in blue. Radial profile averaged over the whole direction (shown in figure 2) is indicated 
by gray line in each panel. 
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Table 1. Results of elliptical fitting in figure 5. 
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inner ring (Biller et al. 2014; Reggiani et al. 2014). Another intriguing feature is that the 
location of deep minimum in the southern part of the inner power-law region seems to coincide 
with the point source at A = 7 mm revealed with J-VLA, as indicated in figure 5b'. Osorio 
et al. (2014) speculate that it might originate from the circumplanetary disk associated with 
a protoplanet inside the gap, which is also identified in thermal emission at A = 7 mm. The 
relation between these features and the possible protoplanet in TZ will be discussed further in 
§4.2. 

To quantitatively describe the shape of TZ, elliptical fitting is made to the contours of 40 
and 80 mJy arcsec -2 in modulated PI; the results are presented in table 1. These ellipses well 
delineate the edges of TZ in the original PI map, as shown in figure 5c. The major axes of the 
ellipses nearly agree with the disk major axis (PA= 5°), and the ratio of minor to major axes 
is ~ 0.8, which is smaller than the case when an infinitesimally thin disk is observed from the 
inclination angle of 13°. These are consequences of the asymmetry that the PI near the minor 
axis of the disk is brighter than near the major axis (see also figures 1 and 3). The centers of 
the ellipses agree with each other but are shifted ~ 07047 northwest from the stellar position, 
caused by the east-west asymmetry along the minor axis in the outer power-law region. 

4. Discussion 

In the following subsections, we discuss observed features of the disk from larger to 
smaller size scales in order; a double power-law distribution seen in the radial profile of PI is 
discussed in §4.1, possible presence of a protoplanet in the transition zone (TZ) is examined in 
§4.2, and asymmetries on smaller azimuthal scales are mentioned in §4.3. It is impossible to 
derive the realistic internal structure of the disk solely from the PI image in 77-band. Existence 
of gap between 40 — 70 AU in radius has been inferred by recent J-VLA observations (Osorio et 
al. 2014), which provided us another important information about the disk internal structure. 
The emission at A = 7 mm, however, is mainly arising from dust grain with a larger size, and we 
cannot make any combined analysis of these two images unless we put further assumptions on 
dust size distribution as well as on size-dependence of dust spatial distribution. We therefore 
decide in this paper to interpret 77-band image with a simple disk model whose surface density 

a Offset from the stellar position 
b Only in 40° < PA <315°. 
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Fig. 5. (a) The image in r < 1"1 of “modulated PI”, which is the PI multiplied by the dimensionless 
parameter r| gAU , where r 2 gAU is the radius normalized by 29 AU. White cross and gray circle shows the 
stellar position and the regions affected by the mask, respectively, (b) The image of PI multiplied by 
r 29AU a ft er the convolution with a Gaussian function whose FWHM is 8 pixels, 0"072. The contours of 
10 and 20 mJy arcsec -2 are shown in white broken lines and white solid lines, and the levels of 40 and 
80 mJy arcsec -2 are shown in black solid lines, (b’) The results of elliptical fittings to the contours of 40 
mJy arcsec -2 (blue) and 80 mJy arcsec -2 except —45° < PA < 40° (red). The centers of the ellipses are 
indicated by blue and red crosses, though the latter is hard to see because they almost perfectly agree 
with each other (see Table 1). The pixels used in the fitting are shown in black crosses, (c) The PI image, 
same as figure la, superposed on the results of the elliptical fittings. 
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and temperature distributions (£ and T) have a power-law and the vertical density distribution 
is determined by hydrostatic equilibrium (Appendix 1). Although our disk model is simple, it 
is useful for analyzing the location of scattering surface, the key to interpret the PI image in 
//-band. As explained in the Appendix 1, S essentially parameterizes the volume emissivity of 
disk material while T parameterizes the vertical distribution of disk material. Further analysis 
on disk structure based on a more realistic model should be made after sensitive observations 
of dust emission at millimeter/sub- millimeter wavelengths are made with ALMA in angular 
resolution similar to this study. 

4-1. Disk structure that reproduces a double power-law distribution of Pi’s radial profile 

An azimuthally-averaged radial profile of the PI in H- band shows a double power-law 
distribution, and the exponents in the both power- law regions are —3 (figure 2). The same 
exponent was also obtained in total intensity at A = 1.1 /iin and 1.6 //m (Grady et ah 2007; 
Fukagawa et al. 2010). A geometrically thin but optically thick disk model (e.g., Whitney 
& Hartmann 1992) has been frequently referred to as an analytical explanation for an re- 
distribution (e.g., Wisniewski et al. 2008), but this is not the only explanation for such radial 
dependence. Grady et al. (2007) showed that a model disk with a pressure scale height of 
H{r ) oc r +1 065 can reproduce the redistribution of surface brightness at A = 1.1 /an by using 
3D Monte Carlo radiative transfer code (Whitney et al. 2003a; Whitney et al. 2003b; Whitney 
et al. 2004). It was also reported that such a minimally flaring disk (i.e. , the aspect ratio 
H(r)/r is independent of r) 2 with surface density E(r) oc r” 1 can explain not only the radial 
distribution of scattered light but also the visibility amplitudes at A = 1.3 mm (Panic et al. 
2008) and the SED at longer wavelengths (Meeus et al. 2010; Honda et al. 2012). 

In fact, r” 3 -dependence of the disk emission may be a natural consequence by the combi- 
nation of E(r) ocr ” 1 and H[r) oc r. The scattering surface of a disk corresponds to the positions 
where the optical depth measured radially from the star is ~ 1, while the energy flux of incident 
stellar radiation can be scaled by (3/r 2 , where fi is the grazing angle at the scattering surface. It 
can be derived analytically that (3 is nearly proportional to r ” 1 for an axisymmetric disk when 
E(r) ocr ” 1 and H(r ) ocr, accounting for the r” 3 -brightness distribution (Appendix 2). This 
means that the double power-law distribution (figures 2 and 5) is realized if fi is proportional 
to r ” 1 in each power-law region but its value in the inner region is significantly smaller than 
that in the outer region. It can happen when the height of the scattering surface normalized 
by the radius in the inner region is smaller than that in the outer region. 

A simple example for such a case is that the surface density distribution of dust particles 
responsible for the scattering in near infrared, Ed(r), is not a smooth power- law but has an 
irregularity. To roughly estimate the amount of irregularity in Ed(r), we compare the observed 
PI distribution with a simple model calculation whose details are described in Appendix 1. As 

2 Such a case has been referred to as a “flat disk” in some studies (e.g., Chiang & Goldreich 1997). 
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discussed later, Wagner et al. (2015) also carry out more elaborate modeling, but our simple 
model is helpful for analyzing the location of scattering surface. £ d (r) of the model disk is 
expressed by 

£ d (r) = / d (r)£ out f — ) \ (3) 

V T , ou t / 

where 

[ fin (r<n n ), 

fd(r) — | hz(r) (r- m <r < r out ), (4) 

[ 1 (r out <r). 

Here, / in s (0 < / in s < 1) is a depletion factor for the surface density in the inner power-law 

region, (r in ,r out ) are the inner and outer boundaries of TZ and are set to be (50 AU, 85 AU) 
based on the fitting results expressed by equations (1) and (2), and £ out is the surface density 
at r out . The surface brightness distribution normalized by r^u does not show any gap but 
smooth rise in TZ (figure 2c), and this cannot be reproduce by a model disk with a dust devoid 
zone. To smoothly connect £ d (r) in TZ as well as at its boundaries, /tzM in equation (4) is 
defined as a third order polynomial that simultaneously satisfies the following equations: 

/tz(Uii) = /in S , 

/Tz(^out) =1, /p-x 

^/tz(Hii) =0, 

J^/TZ^out) =0. 

We also assume that T(r) oc r -1 to set the pressure scale height H{r ) oc r +1 (Appendix 1; see 
also Cliiang & Goldreich 1997), and the disk is isothermal and in hydrostatic equilibrium along 
the vertical axis. The multiple scattering effect is ignored in the radiative transfer calculation 
for simplicity. Figure 6a shows the results of the model calculations. The r~ 3 -dependence is 
reproduced in both r < 50 AU and r > 85 AU regions, but a larger proportionality constant for 
£ d (r) in the outer region makes the location of the scattering surface higher from the mid-plane, 
resulting in larger /3 and brighter PI. Judging from figure 6a, / in s = 1CU 0 ' 8 (~ 0.16) seems to 
give the best fit to the observed contrast between the two power- law regions. £ d (r) given by 
equations (4) and (5) with / in s = 10“ 0,8 is presented in figure 6b. There seems a gap in r f»40 — 70 
AU with the minimum at r & 55 AU as suggested by previous studies (Quanz et al. 2013; Osorio 
et al. 2014), though the profile shape is qualitatively different from a planetary gap having a 
characteristic size scale of the width (e.g., Zhu et al. 2011). Despite many simplifications, the 
estimate by our model is roughly consistent with the results obtained by Wagner et al. (2015); 
/ in s ps 0.3 in the model disk with H(r) oc r +L3 . 

However, the above depletion factor for the inner power-law region should be regarded 
as the lower limit, i.e. , should be > 0.16. This is because irregularity in radial temperature 
distribution can also produce the observed double power-law distribution of PI, even without 
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Fig. 6. Models calculation that reproduce the averaged radial profile ( black crosses ) shown in figure 2. 
Surface density distribution £d(r) and scale height of the model disk H(r ) is proportional to r -1 and 
r +1 , respectively, in r = 29 — 50 AU and r = 85 — 160 AU. £ out and T out at r out = 85 AU are set to be 
9 g cm -2 and 53 K, but these choices are not important in the estimate for / in s or / ir ^ (see Appendix 
for more detail), (a) The proportionality constant for Ed(r) is changed, as described in the equation (3). 
The lines in green, red and blue represent the results of model calculations for /A" = 10~°' 7 , 10~° 8 and 
10“°.9, respectively, (b) £d(r) given by equations (3) and (4) when / in s = lO -08 , normalized by the surface 
density at r = 29 AU. (c) The proportionality constant for T(r) is changed, as described in the equations 
(3) and (4) but replace each E with T. The lines in green, red and blue represent the results of model 
calculations for f ir J = lO -0 07 , 10 -0 08 and 10~ 0 09 , respectively, (d) T(r ) when / it T = lO -0 09 , normalized 
by the temperature at r = 29 AU. 

any irregularity in E^r) like equations (4) and (5); if the temperature proportionality constant 
is slightly higher in the outer power-law region, it makes the disk aspect ratio h = H(r)/r and 
f3 in this region also larger. Our simple model similar to the above one but each £ in equations 
(3) - (5) is replaced by T reproduces the observed PI profile when / in 7 = lO” 0 09 — 10”°'° 8 (~ 
0.81 — 0.83), corresponding to the change in h by a factor of 0.90 — 0.91 (see Appendix 2 for 
further details). This is again consistent with the model calculations by Wagner et al. (2015) 
who show that a scale height change of 0.86 in the inner power-law region can reproduce the 
observed PI image in //-band 3 . There is a degeneracy between / in s and / iu T in models that 
reproduce the PI image, as pointed out by Wagner et al. (2015). The radial intensity profile at 

3 Another important conclusion obtained by Wagner et al. (2015) is that the apparent gap in PI at 40 < r < 70 
AU cannot be reproduced solely by changing the scale height of the inner wall at r « 25 AU. 
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A = 7 mm is reproduced by a disk model in which there is no significant discrepancy in dust 
surface density between r = 30 — 40 AU and r > 70 AU (Osorio et al. 2014), which might favor 
a larger / in s . Further imaging studies, especially at millimeter and submillimeter wavelengths, 
will still be necessary to clearly disentangle this degeneracy. 

4-2. Possible presence of a protoplanet in the transition zone (TZ) 

The change of /3 described in §4.1 may be caused by interaction between a protoplanet 
and the disk. Jang-Condell & Turner (2012), for example, examine the detail disk structure 
under the influence of a forming planet and simulate the images at various wavelengths. They 
obtain a radial distribution of near infrared scattered light qualitatively similar to our observed 
results. Such a double power-law distribution is produced because the protoplanet creates a 
gap in the disk initially having flared surface, resulting in larger /3 on the scattering surface 
at outer radii. Furthermore, cooling within the trough and heating on the far wall alter the 
vertical density profile of the gap region in such a way that the difference of (3 is enhanced even 
more. In their simulations, the brightness contrast between the inside and outside of the gap 
is much larger (a factor of ~ 10 2 in normalized intensity even in the lowest contrast case) than 
our observational results, suggesting that the mass of a protoplanet in TZ is so low that it only 
creates a shallower gap than the cases of Jang-Condell & Turner (2012); they deal with the 
cases of planet masses much higher than the critical mass for a gap opening proposed by Crida 
et al. (2006). 

More recently, Osorio et al. (2014) obtain the image of dust thermal emission at A = 7 mm 
with J-VLA. They detect ring-like emission delineating the bright rim at r « 072 revealed by 
the H - band PI image with VLT/NACO (Quanz et al. 2013). The azimuthally-averaged radial 
intensity profile at A = 7 mm gives a hint that there is an annular gap in r = 40 — 70 AU from the 
star that coincides with the gap seen in the H - band PI images. Based on these results, Osorio 
et al. (2014) carry out disk modeling in which there is a dust devoid zone in r = 40 — 70 AU, and 
they successfully reproduce the radial intensity profile at A = 7 mm. Meanwhile, the scattered 
light in 77-band is detected even in the gap (figures 1 & 2), indicating that there is some amount 
of small dust particles in this zone. These two observational results are compatible if larger 
dust particles responsible for thermal emission at millimeter wavelengths are more deficient in 
the gap region. In fact, previous studies on the motion of dust particles in a gas disk under the 
influence of a planet show that particles with a modest size (i.e. , millimeter and centimeter in 
size) tend to be removed more quickly compared to the smaller particles during the formation 
of a planetary gap near the orbit (Paardekooper & Mellema 2006; Lyra et al. 2009; Muto & 
Inutsuka 2009). 

The J-VLA image at A = 7 mm also reveals a faint point source, 0734 at PA=175° from 
the central star. Osorio et al. (2014) discuss the possibility that it arises from thermal dust 
radiation of a circumplanetary disk associated with the protoplanet sculpting the gap. In the 
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image of modulated PI in H- band, on the other hand, there is a deep minimum in the inner 
power-law region whose location seems to coincide with the point source at A = 7 mm (figure 
5b'; see also figure 3a). This can be regarded as another sign of the protoplanet, because the 
disk scattering surface can be lowered by the protoplanet’s gravity, making an observable dip 
in near infrared scattered light (Jang-Condell 2009). 

To constrain the planet mass, it is important to precisely measure the depth and width 
of the gap and to compare them with a latest theoretical model of gap structure (e.g., Duffell 
& MacFadyen 2013; Kanagawa et ah 2015a). An analytic formula among the gap depth, planet 
mass, disk aspect ratio and viscosity is presented in a separate paper (Kanagawa et ah 2015b), 
in which the formula is applied to the gap depth inferred from J-VLA observations of HD 
169142 by Osorio et ah (2014). Imaging studies of both dust and gas emission with the ALMA 
will also be important in near future. 

4-3. Asymmetric features in the outer power-law regions 

As mentioned above, there are remarkable asymmetries both in the inner and outer 
power-law regions. Significant asymmetric variation in the inner power-law region (figures 3 
and 5) could be due to the perturbation by a protoplanet inside the inner ring (Biller et ah 
2014; Reggiani et ah 2014). The deep minimum of PI in H - band in the inner power-law region 
coincides with the point source at A = 7 mm, which can be regarded as another sign of a 
protoplanet surrounding the circumplanetary disk. The azimuthal distributions at the outer 
radii tend to have the maxima at PAf« 110°, along the disk minor axis east to the star. 

Besides these features, there are small-scale asymmetries in the outer power-law region 
(see figure 5), though the amount of variation is smaller than in the inner power- law region. 
There are two possibilities for the origin of these small-scale asymmetries. One is corrugation of 
the scattering surface in the outer region, causing a small-scale variation in grazing angle (/3). 
Such corrugation may be due to local fluctuations in disk surface density that may affect the 
self- gravity of the disk and hydrostatic balance in its vertical direction (Muto 2011). The other 
possibility is a shadowing effect of the structure in the inner power-law regions. As described 
in §4.1, the radial intensity profile of near infrared scattered light as well as the broadband 
SED are basically explained by a minimally flaring disk (see also Grady et al. 2007; Panic et al. 
2008; Meeus et al. 2010; Honda et al. 2012; Wagner et ah 2015). If a part of the inner structure 
of such a disk is slightly puffed up, it may cast a shadow on the outer region observable in the 
near infrared image. Radial profiles along the major and minor axes (figure 4) actually show 
that the surface brightness in the inner (r < 50 AU) and the outer (r > 80 AU) regions are 
anti-correlated with each other. Similar tendency is also seen in all the pairs of radial profile 
180° apart in PA from each other (figure 7). These two possibilities may be disentangled by 
search for time variation on a shorter timescale in scattered-light images; if the asymmetries in 
the outer disk are caused by shadowing of the inner region, asymmetric patterns in the outer 
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Fig. 7. Radial profiles of the PI along the pairs separated by 180° in PA. (a) PA = 5° ± 5° in red and 
PA = 185° ± 5° in blue. Same as figure 4a. (b) PA = 95° ± 5° in red and PA = 275° ± 5° in blue. Same 
as figure 4b. (c) PA = 35° ± 5° in red and PA = 215° ± 5° in blue, (d) PA = 65° ± 5° in red and PA = 
245° ± 5° in blue, (e) PA = 125° ± 5° in red and PA = 305° ± 5° in blue, (f) PA = 155° ± 5° in red and 
PA = 335° ± 5° in blue. Radial profile averaged over the whole direction (shown in figure 2) is indicated 
by gray broken line in each panel. 

disk may vary on a short timescale that synchronizes with the change of disk structure in the 
inner power-law region (see also Wagner et ah 2015 for discussion on the relationship between 
the time variation seen in SED and structural change at much inner radii as small as r ~ 0"3). 

The PI in the outer region tends to be higher along the minor axis than the major axis 
(figure 4). This is completely opposite to other disks showing brighter PI along the minor 
axis (Muto et al. 2012; Kusakabe et ah 2012), which is naturally explained by the highest 
polarization degree in 90° scattering. The observed features of the disk around HD 169142 are 
rather difficult to interpret in terms of simple Rayleigh scattering by an axisymmetric disk in 
nearly face-on geometry, ft is more likely that the most of these asymmetric features reflect 
real disk structure on smaller scales, as discussed above. 

5. Conclusions 

We carry out coronagraphic observations of the circumstellar disk around HD 169142 in 
H- band polarized intensity (PI) with Subaru/HiCIAO. The emission scattered by dust particles 
at the disk surface in 0"2 < r < H'2, corresponding to 29 < r < 174 AU at the distance to HD 
169142, is successfully detected. Our conclusions are summarized as follows: 
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1. The overall structure of the PI is nearly axisymmetric. The azimuthally- averaged radial 
profile of the PI shows a double power- law distribution, in which the Pis in r = 29 — 52 AU 
and r = 81.2 — 145 AU respectively show r” 3 -dependence. These two power- law regions 
are connected smoothly with a transition zone (TZ). The apparent gap in r = 40 — 70 AU, 
previously reported by Quanz et ah (2013) and Osorio et al. (2014), corresponds to the 
outer part of the inner power-law region and TZ. 

2. The PI image also contains asymmetric features. The variation along azimuthal direction 
in each annulus is systematically larger at inner radii, and the amount of variation in the 
inner power-law region is greater than 4. This could be due to the perturbation by a 
protoplanet inside the inner ring (Biller et al. 2014; Reggiani et al. 2014). The azimuthal 
distributions at the outer radii tend to have the maxima at PA^ 110°, implying that the 
eastern part is the near side of the disk and is brighter than the western part because of 
the forward scattering. 

3. The radial profile of r _3 -dependence is realized by the combination of surface density 
E(r) oc r^ 1 and the pressure scale height 77(r) oc r. The observed double power-law distri- 
bution of the PI can be reproduced when the proportionality constant of surface density 
in the inner power-law region is smaller than that in the outer power-law region. Simple 
model calculation shows that the observed PI profile matches the case of / in L = 0.16, where 
/ in s is the depletion factor for the proportionality constant of surface density in the inner 
power-law region. However, this should be regarded as a lower limit because irregularity 
in temperature can also produce a double power-law distribution. Even with no irreg- 
ularity in surface density, the observed PI prohle is reproduced when / in T = 0.81 — 0.83, 
where / in r is the proportionality constant of temperature in the inner power-law region 
relative to the outer region. These results are consistent with more elaborate modeling by 
Wagner et al. (2015). The degeneracy between / in s and / in T in models that reproduce the 
PI image should be disentangled by further imaging studies, especially at millimeter and 
submillimeter wavelengths. 

4. The modeling based on the image at A = 7 mm (Osorio et al. 2014) requires a dust devoid 
zone between r = 40 — 70 AU, while the scattered light in 77-band is significantly detected 
even in this gap. These observations suggest that larger dust particles responsible for 
thermal emission at millimeter wavelengths may be more deficient in the gap, as predicted 
by previous theoretical studies. 

5. The PI in the inner power-law region shows a deep minimum whose location seems to 
coincide with the point source at A = 7 mm. This can be regarded as another sign of a 
protoplanet, because the scattering surface of the disk can be lowered by the protoplanet’s 
gravity, making an observable dip in near infrared scattered light (Jang-Condell 2009). 

6. The obtained PI image also shows small scale asymmetries in the outer power-law region. 
Possible origins for these asymmetries include corrugation of the scattering surface in the 
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outer region, and shadowing effect by a puffed up structure in the inner power-law region. 
These two may be disentangled by search for time variation on a shorter timescale in the 
scattered-light image. 
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Appendix 1. Description of the model calculations 


The model is based on a disk where the surface density and temperature distributions 
are given by power-law forms: 
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where E out and T ou t are the surface density and temperature at r out . In the model, the sound 
speed c s is approximated by 
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for simplicity. 

The disk is assumed to be vertically isothermal and in hydrostatic balance with the 

gravity of the central star. The density profile is then expressed by 
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where p c (r) is the midplane density and H = c s /f2 is the pressure scale height. Assuming the 
Keplerian rotation, it is easy to see that H ocr^ 3-9 ^ 2 is realized from equations (A2) and (A3). 
Using equation (A3), the disk aspect ratio h = H/r scales with 
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The midplane density, on the other hand, is related with the surface density as 
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Since our goal in this study is to assess the location of the scattering surface in a simple setup, 
we do not care about the consistency between the temperature and density profile in terms of 
the energy balance. 
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In radiative transfer calculation, we assume a single and isotropic scattering. 1 The 
equation to be solved is 

= - a w(r, z)B(T*) exp[— r rad (Z) - t(Z)\ , (A7) 

where a is the scattering coefficient, Z is the coordinate along the line of sight, W(r,z) = 
l/4(r 2 + z 2 ) is the geometric dilution factor, r ra d is the optical depth between the central star 
and the location of the disk, and t(Z) is the optical depth along the line of sight. In integrating 
equation (A7), we first calculate r ra d at all the location of the disk in the grid and then integrate 
over line of sight from the observer position to the disk. 

To compare with the observed radial profile of the PI, we modify the surface density or 
temperature distribution as described in §4. The averaged PI is assumed to be proportional to 
the averaged total intensity, which seems reasonable because the total intensity has a similar 
radial profile to PI (Grady et al. 2007; Fukagawa et al. 2010). We treat the unit of calculated 
intensities as an arbitrary one partly because of the uncertainty in polarization efficiency, and 
we fit the shape of the PI radial profile to estimate / in s or fZ. 

In this paper, we only deal with the cases of p — 1 and q= 1 and look at how the observed 
radial profile can be reproduced. The q = 1 case corresponds to the disk whose aspect ratio 
h is independent of r as shown in equation (A5) (Grady et al. 2007; Panic et al. 2008) and 
successfully reproduces the r~ 3 -dependence of the PI profile. The r _3 -dependence can also be 
derived analytically in p = 1 and q — 1 cases, as shown in Appendix 2. Figure 6 shows the 
cases of £ out = 9 g cm -2 , T out = 53 K at r out = 85 AU taken from equilibrium temperature for 
blackbodies directly exposed to the stellar radiation with 9.4 L & (Meeus et al. 2012). Although 
the fiducial grain temperatures are set to the local blackbody temperature, the grains are 
assumed to be highly reflective, with the scattering albedo of uj = 0.9, and their total extinction 
coefficient is set to be x — 100 cm 2 g" 1 (a = yu>). The change of £ out or T out , however, almost 
scales the overall profile unless the disk gets optically thin, and it does not significantly affect 
the estimate for / in s or / in T . 

The position of the scattering surface depends on two factors: (i) the volume emissivity 
of disk material determined by its density and emissivity, and (ii) the vertical distribution of 
disk material. The model described here essentially parameterizes (i) with £ and (ii) with T, 
respectively. The change of £ can be interpreted not only by the change of material density 
but also by the spatial variation of optical properties of dust particles. The change of T, on the 
other hand, can be related not only to the temperature, or the pressure scale height H = c s /f2, 
but also to the degree of dust settlement toward the disk mid-plane. Since it is impossible 
to disentangle all these factors solely from the PI image in H- band, we do not examine the 
detailed disk structure based on a realistic disk model in this study. Further detailed analysis on 

4 Thermal emission is included in the code as well, but it is not significant as long as we consider the intensity 
in near infrared beyond 10 AU. 
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disk structure should be made after sensitive observations of dust emission at millimeter /sub- 


millimeter wavelengths are carried out with ALMA in a similar angular resolution. 

Appendix 2. Analytic consideration for the cases of p = 1 and q — 1 

The results of the model calculation show that the cases of p = 1 and q = 1 successfully 
reproduce the r~ 3 -dependence of surface brightness. In fact, it is possible to predict the radial 
profile of the surface brightness analytically in these cases as follows. 

The approximate solution to equation (A7) is given by 


where /3 is the grazing angle of the stellar light at the disk surface. Here, the disk scattering 
surface is determined by the height z s at which the optical depth from the central star becomes 
unity. The optical depth r ra( j can be calculated by 
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where 6 is given by 


tan t) = 
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Substituting the density profile given by equation (A4), we have 
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Here, we have assumed that the central star is at the origin and used the fact that z' = r' tan# 
along the line connecting between the location of the disk considered and the stellar position. 
If we assume the disk with the constant h, the exponential term is constant in the integration 
and therefore 
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Finally, if the midplane density profile is given by 
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its integration can be calculated as 
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Note that r ra d(r, z) has the form 

Tnd(r,z) = f(r)g(z/r) = f(r)g(u), 
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We can define the disk scattering surface z s (r ) by 

r^d{r,z s ) = t 0 , (A18) 

where To is the value of the order unity. Then, the location of the scattering surface is given by 
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where we have used h <C 1 and 
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In the regions of r ^ r in , z s /r in equation (A19) only weakly depends on the radius, and the 
scattering surface always resides a few times higher than the disk scale height. Unlike a flared 
disk (e.g., Chiang & Goldreich 1997), the aspect ratio of the scattering surface, z s /r, is almost 
constant (i.e. , the scattering surface is nearly parallel to the radial direction from the star; see 
equation (A19) in Appendix 2) when p = 1 and q = 1. In such a case, the surface brightness of 
the resultant scattering light is changed rather sensitively by the change of £. 

The grazing angle can be calculated by 
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where ' denotes the derivative with respect to the argument. Substituting the expressions of / 
and g functions and use h • Cl, we obtain 
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In r ri n , the last part in a bracket is a weak function of r. We therefore conclude that the 
grazing angle scales with 

/3 oc r~\ (A25) 

and combining this with equation (A8), the radial profile of the scattered light is given by 

/ oc r -3 . (A26) 

Equation (A24) indicates that f3 is nearly proportional to h when £ distributes smoothly, 
while h is proportional to c s , or T 1 / 2 . The difference in observed proportionality constant of 
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PI(r) between r < 50 AU and r > 85 AU is 0.22, as shown in equations (1) and (2). This 
difference corresponded to 0.47 in temperature if the aspect ratio h remained constant at all 
radii. In our model, however, h in r < 50 AU is smaller than that in r > 85 AU, and the 
stellar radiation incident at the scattering surface at r > 85 AU is brighter than the case of h = 
constant at all radii. This is why the models with / in T closer to the unity (= 0.81 — 0.83) in 
§4.1 reproduces the observed radial profile. 
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